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O : ABSTRACT 
(N . 

, ^ . Calculations are performed for the spectral line profiles and images of as- 

Q , trophysical maser radiation that emerges from isolated spheres and thin disks 

i viewed edge-on. In contrast to previous investigations where various approxima- 

^ ! tions are made, the full equations are solved here for the frequency-dependent 

^ ! radiative transport that includes the thermal motion of the molecules. The spec- 

Tj- ! tral line profiles for spheres and disks are found to rebroaden to the full thermal 

^ ! Doppler breadth with increasing saturation in essentially the same way as is well 

O ! known to occur for a linear maser. The variation with frequency in the apparent 

Q ! angular sizes of masing spheres and thin disks is found to be negligible at fre- 

r-| ! quencies within the spectral line where the flux is significant. Calculations also 

^! are performed for spherical and disk masers that are not isolated, but for which 

P ! the seed radiation for the masers is incident from one side as would occur when 

i a strong continuum source is on the far side of the masers. Again, the spectral 

i line profiles are found to rebroaden to the full thermal breadths with increasing 

.!-H i saturation and there are no significant variations in the apparent angular sizes 

^ I with frequency. However, the full rebroadening does occur at somewhat higher 

5t ! saturation and the variation of the apparent angular sizes as a function of the 
degree of saturation is quite different from that of the isolated masers. Spheres 
and disks have served as idealized geometries with which to examine possible 
deviations from the linear approximation for astrophysical masers. 

Subject headings: ISM : molecules — line profiles — masers — radiative transfer — 
radio lines : ISM 



1. INTRODUCTION 



Interpretations of astrophysical masers which involve the radiative transfer of the maser 
radiation are based almost entirely on the "linear maser" (i.e., the one dimensional) ap- 
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proximation. In this approximation, the paths of the rays of radiation within the maser are 
parallel, straight lines. This approximation is motivated by the recognition that the longest 
maser rays tend to dominate as a result of the exponential gain in the intensity with path 
length in the unsaturated regime. From an observational viewpoint, that the rays are nearly 
parallel is indicated by evidence that the angles into which the emergent maser radiation is 
beamed are quite small — as small as 10~^ ster in some cases (Nedoluha & Watson 1991). 

As an alternative to the extreme approximation that the masers are purely linear, the 
transfer of radiation in isolated spherical masers has been treated in a few investigations. 
Although it is generally recognized that astrophysical masers are unlikely to be ideal spheres, 
spheres nevertheless serve as a convenient geometry with which to assess possible deviations 
from the linear maser approximation. The initial studies (Litvak 1971; Goldreich & Keeley 
1972) demonstrated that the maser radiation becomes tightly beamed along the radii so 
that the rays are again nearly parallel. Consequently, with increasing saturation the apparent 
sizes of these masers as seen by a distant observer become much smaller than the geometrical 
sizes. In these studies, the maser radiation was treated as containing only a single frequency 
and the thermal (and any other) velocities of the molecules were ignored. Efforts to treat 
the full spectrum of the maser radiation and the thermal motions of the molecules in a 
spherical maser seems to have begun with Bettwciscr & Kegel (1974). The key aspect of the 
radiative transfer for masers that are not completely linear was recognized — the frequency 
dependence of the amplification of a ray varies with its direction of propagation. For spheres, 
the amplification depends upon the angle between the ray and the radius vector. Bettweiser 
(1976) concluded that the observable properties of spherical masers, such as the spectral 
line breadths and the apparent sizes as a function of frequency, may be sensitive to the 
details of the treatment of the frequency-dependent radiative transport. However, with 
limited computer capabilities, definite conclusions could not be reached in that investigation. 
To simplify the computations somewhat, Emmcring & Watson (1994; hereafter EW) then 
developed a series expansion method that is based on the knowledge that the radiation is 
beamed into small angles about the radii. Despite the complications in the detailed treatment 
of the anisotropic, frequency dependent radiative transport, on the issues of main interest 
EW found (a) that the spectral line narrows and rebroadens as a function of the degree of 
saturation in essentially the same way as for a linear maser, and (b) that the variation in 
the apparent angular size as a function of frequency is negligible at frequencies within the 
spectral line where the flux is appreciable. Speculation to the contrary on these issues has 
nevertheless persisted — even including assertions of a specific behavior that the line profiles 
from disks and spheres rebroaden only to l/\/2 and l/V^ of the thermal Doppler breadth, 
respectively, instead of rebroadening to the full Doppler breadth. The calculation of EW 
does employ an involved expansion proceedure and some reasoning is required to recognize 
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that the results are accurate. The reasoning may not be completely clear to all readers. 
To remove these challenges to understanding that the key conclusions are valid, we have 
now computed the spectra and images by solving the full equations without approximations. 
In addition to spherical masers, masers in the shape of thin disks and viewed edge-on also 
are treated. Disk masers are, perhaps, more likely than spherical masers to occur in the 
astronomical environment. The results from the calculations with the full equations are 
described in this Paper. 

The issues under investigation here are relevant for interpreting observational data. The 
observed breadths and shapes of the spectral lines of astrophysical masers can be used to 
infer information about the temperature, the degree of saturation, and the beaming angle of 
masers (e.g., Nedoluha & Watson 1991; Watson, Sarma, & Singleton 2002) if the theoretical 
behavior of the profile is understood. It is striking that the spectral hnebreadths associated 
with some of the brightest astrophysical masers (the H2O maser flares) are much less than 
the Doppler breadths for the temperatures (> 400K) that are expected for the gas in the 
relevant environments. The images of astrophysical masers often are not identical when 
observed at different frequencies within the spectral line (e.g., Gwinn 1994; Moscadelli et al. 
2003). Establishing whether this difference is, at least in part, a natural consequence of the 
radiative transfer is necessary when considering its significance. 

The relevant equations to be solved are given in Section 2. Results of compTitations for 
the spectral line breadth and for the apparent angular size of the masers arc presented in 
Section 3. As in the previous investigations, the focus in the above discussion is on "isolated" 
masers — that is, on masers for which the seed radiation is due either to spontaneous emission 
or to isotropic continuum radiation that is incident from outside the masing region. In 
practice, some masers in the astronomical environment may be close enough to HII regions 
or other strong sources of continuum radiation that the radiation from one side dominates 
as the source of the seed radiation. Calculations are performed in Section 4 for this 
well, and compared with the spectral line breadths and apparent angular sizes obtained for 
isolated spherical and disk masers. 

2. BASIC EQUATIONS 

For simplicity, we consider a masing transition between two energy levels which have 
statistical weights of one. Let riill^ ,lf')d'r^dlf be the number of atoms of level i with 
velocities in the range {if, if + dlf) and located in the volume element dl^ at position 1^ . 
The rate equations for the populations of the upper u and lower / levels can then be written 
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as 

= Xu{v) - Tuui'T^, if) - lf)BJk{'r^, if) - Auuil^, if) (1) 

and 

= Xi(v) - TniiT', If) + n_( r", ir)BJk(T, if) + r", if) (2) 

Here, the Aj are the pump rates which always are assumed to be proportional to the 
MaxweUian velocity distributions, F is the loss from the states due to all processes other 
than the maser transition, A and B are the Einstein coefficients which are related by 
B — c^A/2hu^ since the statistical weights of the upper and lower states are assumed to be 
equal, 

K{v) - Xi{v) 1 
— Uu — m — ^ ^ , , ,^ (6) 

r i + 2BJk{i^,if)/r ^ ' 

when an unimportant term in {A/T) is ignored for simplicity and 

Jk{y,lf)= J I[y,k,iy^iyo{l + k-lf/c)]dk/in (4) 

Here, I{l^,k,i/) is the intensity of a ray of maser radiation propagating in the direction k 
at a location with frequency z/. The increase in / when it traverses a differential path 
length ds is given by the radiative transfer equation (EW) 

dl{l^,k,i')/ds = (huo/ATi) J G{1^, if ,k, i')d[i' — uq{1 + k ■ if /c)]dlf (5) 

where G{1^, if, k, u) = n_(7^, lf)BI(T^, k, v) + Auui/f , if). It is convenient to express the 
foregoing equations in terms of normalized quantities I, Jk , and k so that they become 

dl{l^, k, u)/dr — k{1^, k, u)I{l^, k, v) + ?7exp(— i?^) (6) 

where dr is the differential optical depth parameter and the normalized maser opacity 

TT y 1 + v!) 

is given in terms of the integral over the normalized velocities l6 of the masing molecules 
and 

Jk{r,w,9') = — J sin 9d9 J I {r, 9, u = w cos 9 cos 9' + w sin 9 sin 9' cos (j))d4> (8) 

which expresses the rate of stimulated transitions for molecules with a normalized velocity 
of magnitude w in the direction given by the angle 9' measured relative to the radius vector. 
In writing these expressions for k and Jk, we have recognized that the calculations here 
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will only be applied to geometries (sphere and disk) other than the linear maser where 
the intensity only depends on the radial distance from the center and the angle 6 between 
this radius and k. The dimensionless I{l^,k,u) and Jk{l^,ll}) in the above equations and 
throughout the remainder of this Paper arc obtained by dividing the actual intensity by 
a characteristic "saturation intensity". In the above equations, the saturation intensity is 
Is = hv^T / (^A. The normalized frequency D is related to the actual frequency v of the 
radiation by z/ = {c/vj:,)^^ — foj/fo where vq is the rest frequency of the transition and 
I'd is the thermal breadth that appears in the Maxwellian distribution of the velocities, 
and is reflected in the pumping rates, e.g. \u{v) — Au(0) exp[— (i>/i>D)^]). The normalized 
molecular velocity is then id — it /vd- In the integrand in equation (8), the frequency v also 
depends upon the anglular difference between the azimuthal angles of k and measured 
about as the axis. The optical depth parameter dr = n^ds where = hcvjjB\_{0)/Ar. 
For simplicity, negligible terms involving the ratio (A/T) have been omitted (see EW). The 
contribution of spontaneous emission in equation (1) is representated by ?7exp(— z/^) where 
77 = (/l/r)[A+(0)/A_(0)] and A+(0),A_(0) arc the sum and difference of the pump rates 
into the upper and lower of the masing states at zero velocity, respectively. The for the 
linear and disk masers depend upon the assumed beaming angles and it is convenient to use 
other expressions for the saturation intensities in the calculations for these geometries. The 
geometry of the calculation is explained in more detail in Appendix B. 



2.1. Linear masers 

For a linear maser, the foregoing equations are greatly simplified and, in the absence of 
spontaneous emission, the radiative transfer equation adopts the familiar form (e.g., Goldre- 
ich & Keeley 1972) 

dl{s, ±ko, V)/dT = ± }y \' 5; ^ + 77 eM-^) (9) 

1 + I{s, +ko, u) + /(s, -fco, I') 

where /(s, +/co) i^) is the intensity that is propagating in the direction of increasing s measured 
along the axis of the maser and I{s, —ko,u) propagates in the opposite direction. This 
dimensionless intensity is obtained by dividing the actual intensity by the saturation intensity 
for a linear maser Is = AnhuQT /c^AAQ. The meaning of the spontaneous emission parameter 
also is modified to incorporate the beaming angle to give r] = (AAfi/47rr)[A+(0)/A_(0)] in 
equation (9). Equation (9) follows from the foregoing general expressions when the integral 
over solid angles in equation (4) is performed by treating the intensities as constant within 
a solid angle equal to the beaming angle AQ of the maser radiation. This is the standard 
relationship between and / in the linear maser idealization. 
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2.2. Spherical masers 

For the spherical maser, one of the integrations for K{r, k, v) can be performed immedi- 
ately to give 

oo TV 



K(r,e,u) = - zdz ^ —di/j (10) 

where w'^ — i^ + z'^ and 9' — arccos[(2;cos'0sin^-|-i/cos^)/u»] (see EW, though note that the 
variables in EW are not exactly the same as here). Physically, ijj is the angle of projection 
of Uf on the plane perpendicular to k . The remaining integrations are performed most 
conveniently in these new variables. Equation (10) is then used together with equations (6) 
and (8) to perform the numerical computations for the spherical maser. 



2.3. Disk masers 

To perform the integral in equation (8) over the angle that is perpendicular to the plane 
of the thin disk and to relate the result to a physical intensity which must have a finite 
angular extent, an idealization is made that is analogous to that for the linear maser. The 
intensity is treated as constant within the beaming angle Aa measured perpendicular to the 
plane of the disk. In the two dimensional geometry that is appropriate for an infintesimally 
thin disk, the integral for then becomes 

2-lT 

Jk{r,w,e') = — / 7(r,6',i? = w cos 6^ cos 6^' + ^ sin 61 sin 6'')(i6' (11) 



when the intensity / is expressed in units of a saturation intensity that is appropriate for the 
two dimensional geometry Ig = 2hvQV / c? AAq. Incorporating Aa into the definition of the 
saturation intensity is analogous to incorporating Af2 into the definition of the saturation 
intensity for the linear maser. With k restricted to be in the plane of the thin disk, two of 
the integrations in equation (7) for fi;(r, k, v) can be performed immediately to give 

oo 

K{r,e,u) = - [ dzeM-y^% — — — + - — 777] (12) 

in which w"^ = 1^ + z"^ and = ^ arccos(z//w) to facilitate the final integration. Again, 9 
is the angle between the radius vector and the direction of propagation, and 9' is the angle 
between the radius vector and the molecular velocity. The computations for disk masers are 
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then performed with these two equations and with equation (6). The term for spontaneous 
emission must also be interpreted as r/ =(AQ;/2)(A/r)[A+(0)/A_(0)] in order to incorporate 
the integration over the angle perpendicular to the plane. That is, equation (6) for the 3D 
geometry is effectively multiplied by Aq;/2 to express it in units that are convenient for a 
disk. 



3. RESULTS OF COMPUTATIONS 

"Saturation" reflects the importance of maser radiation for altering the population dif- 
ference between the masing states, and hence the effect of the radiation upon k. In general, 
this depends upon the velocity Uf of the molecule, and upon the direction and frequency of 
the radiation. As the benchmark for the degree of saturation of a maser, we adopt "R/F" = 
Jk{Re, 0, 0)/2 for disk and spherical masers as in EW. This quantity is the rate for stimulated 
emission (divided by the loss rate F) for molecules with zero velocity at the edge of the masing 
region where the radius is Rg . In practice, our computations for Jk{Re, 0, 0)/2 differ negli- 
gibly from {K{Re, 0, 0)""^ — l}/2 which is a possible alternative, and perhaps more precise, 
measure of the degree of saturation. The analogous expression for the degree of saturation at 
the end of a linear maser of length L is R/F = [7(L, +ko, 0)+I{L, -ko, 0)]/2 ~ 7(L, +ko, 0) /2 
where ko is in the direction of the emerging radiation. 

The equations are solved by iteration. In each step, a new Jk is chosen based on the Jk 
that are obtained in the two most recent steps. This continues until J and k at all w, and 
angles in the computational grid change by less than a specified fractional amount between 
successive iteration steps. We have, of course, verifed that the results that we present are 
insensitive to this fractional difference required for convergence, as well as to the number 
of grid points in w, r, and in the angles 6 and 6' at which J and k are computed. In 
the computation for the disks, J and k, typically are computed at 33 values of each of the 
variables w, z/, r and the angles 9 and 9'. For a sphere, we are unable to allow 33 values 
for all variables if the computations are to be performed in a "reasonable" time of a few 
days. Some variable must be constrained to 17 values, but we verify that the accuracy of 
the computations is not significantly compromised. These values (the "grid" ) are uniformly 
spaced, though the ranges of angles that are included are reduced to reflect the beaming 
of the maser radiation. The ranges of angles depend upon the location in radius at which 
J and K are being computed, as well as upon the degree of saturation of the model being 
considered. Additional details about the geometry and about the computations are given in 
Appendix B and Appendix C. 

Two quantities must be specified at the outset of a computation to solve the foregoing 
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equations for any of the three geometries — 

a) the size of the maser expressed in terms of the number of unsaturated optical depths, 
that is, KgRe for the disk and spherical geometries and KqL for the linear maser (where L is 
the length of the maser) and 

b) the seed radiation parameter — either the external continuum intensity /q that is 
incident onto the outer surface of the masing region or the spontaneous emission term rj. 
We find that the results of computations performed with a specific external intensity Iq and 
zero spontaneous emission differ negligibly from those performed with zero external intensity 
and with the spontaneous emission term rj set equal to /q. In fact, they are indistinguishable 
except at the very highest R/T for which we present the results in the Figures. Hence, we 
will give explicit results only for 77 = and various values of Jo, but with the understanding 
that the results are applicable when spontaneous emission is the dominant source for the 
seed radiation and rj takes on these same numerical values. 

To compare the behavior of the maser properties between the three geometries, the 
natural comparison is between computations for which the background radiation is the same. 
The background radiation is the same when the rates for stimulated emission (with ?7 = 0) are 
equal when no masing gas is present — that is, when R/T is the same for the three geometries 
in the limit where KoRg and kqL = 0. Hence, we compare computations for which Iq is the 
same for the sphere and the disk, and those for a linear maser for which the Iq is one-half of 
that for the sphere and disk. Once Iq is specified, the intensity of the radiation that emerges 
from the maser then depends only upon KoRe or kqL. Instead of designating our results 
by the values of koRc and kqL, we use the degree of saturation R/T which is a measure of 
the emergent intensity. To summarize, the computations will be labelled according to their 
values of R/T (which substitutes for KoR^ or KqL) and Iq. 

What are representative values for /q? Consider the 22 GHz masing transition of the 
water molecule which is observed and imaged extensively in a wide variety of astronomical 
environments. Typical calculations (e.g., Anderson & Watson 1993; Wallin & Watson 1997) 
for the pumping rates yield A+(0)/A_(0) ~ 10 — 100. For r(~ 1 s~^) which is considered 
representative for these masers, then rj = (74/r)[A+(0)/A_(0)] ~ 10~^ to 10"*^. For spherical 
masers, no additional uncertainty enters due to lack of knowledge about the beaming angle. 
Hence, Iq = {A/T){2kTb/hiyo) ~ 2TbA/r ~ 4 x IQ-^ Tb for isotropic background radiation 
with brightness temperature T^. The effective value of Ti, is expected to be close to the 3K 
of the cosmic background. For the hnear and disk geometries, these expressions must be 
multiplied by the beaming angles — by An/47r for the linear geometry and by Aa/2 for 
the disk geometry — to obtain the analogous estimates for Iq and rj for these geometries. 
Beaming angles for astrophysical masers are highly uncertain. Ratios of ten for the length- 
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to-width commonly are discussed (for which AQ/47r ~ 10~^ and Aci;/2 ~ 10~^), though 
the actual beaming angles may easily be smaller by a factor of ten or more. With the 
understanding that the spontaneous emission term can be replaced by an Jq of the same 
magnitude without changing the results, we can then focus our attention on computations 
with Jo < 10^^, at least for applications with the 22 GHz water masers. Analogous estimates 
suggest that this range is representative for other strong masers of methanol and 18 cm 
hydroxyl, though the range may be extended up to ~ 10~^ for the circumstellar SiO masers. 

3.1. Spectra 

The flux of radiation at frequency V from a masing sphere that is seen by a distant 
observer is 

F{v, To) = 2f^geom / I{Re, 0, i") COS 6 siu 6(16 (13) 
Jo 

whereas that from an edge-on disk is 

F{U, To) = Qgeom / ^ I {Re: 0, IP) COS 6d6 (14) 

and that from a linear maser is simply 

F{D, To) = ^geomH^, l^) (15) 

in terms of the intensities at the edges which propagate at angles 6 relative to the radius 
vector (or along in the case of linear masers), and where figcom is the geometrical angular 
size of the maser. Note that the flux F is in units of the saturation intensity for the 
relevant geometry, and that the unsaturated optical depth parameter To = 2«;o-Re or tq = HqL 
depending upon the geometry. 

The behavior of the spectral line breadth of the maser flux as a function of the degree 
of saturation is one of the two main results of this investigation. Representative results are 
given in Figure 1. There the spectral line breadth Vi/2 — obtained by solving -F(??i/2,To) = 
F(0,To)/2 — for the maser flux from disks and spheres is seen to behave in essentially the 
same way as that from linear masers. As a function of the degree of saturation i?/r, the 
breadth first decreases as the maser becomes more intense but is still unsaturated. Then, 
when saturation can no longer be ignored as R/V approaches one, the profile rebroadens and 
ultimately reaches the full thermal Doppler breadth just as occurs for linear masers. The 
rebroadening begins slightly later as a function of R/T for the disk and the sphere. Once it 
begins, however, the broadening is somewhat more rapid than for linear masers. 
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Although the maser profiles only approach the Gaussian that has the breadth of the 
thermal gas velocities when the maser is saturated, the profiles can be quite close to Gaussians 
with smaller breadths in the unsaturated regime. A useful quantitative measure of the 
deviation of the maser profile from a Gaussian is (Watson et al. 2002) 

^ _ I[F{^,ro) - ai^exp{-lP/a2)]'^du 

where Fp = F(0, Tq) is the peak value of the flux which occurs at the center of the line 
profile, 2z/i/2 is the full width at half maximum (FWHM) of this profile, and (01,02) are 
the parameters obtained in fitting F{u, tq) to a Gaussian using the least squares method. 
Although 6 can be small in both the saturated and unsaturated regime, the nature of the 
deviations from Gaussian in the two regimes is different. The computed profile is more 
sharply peaked than a Gaussian (positive kurtosis) in the unsaturated regime, but is less 
peaked than a Gaussian (negative kurtosis) in the saturated regime. The deviation 5 and 
the kurtosis K, where 

K ^ j F{v,TQ)v^dv/a^ j F{v,TQ)dA) (17) 

are shown in Figure 2 for the three geometries. Here, a is the standard deviation of the dis- 
tribution F{i/,tq). Certain observational evidence (Watson et al. 2002; Moscadelli, Menten, 
Walmsley, & Reid 2003) indicates that the 5 and K for actual spectral fines are consistent 
with the predictions in Figure 2 for spectral fines that have not been rebroadened by the 
effects of saturation. 

As for linear masers, tfie linebreadtfis for tfie disks and spfieres tend to increase witfi 
Iq when the maser is not fully saturated as seen in Figure 3. However, the behavior is 
relatively insensitive to Iq. In all cases, the fine breadths are less than the thermal breadth 
and approach the thermal breadth as i^/F exceeds one. 



3.2. Appcirent anguleir sizes 

We calculate the apparent angular size as the region from which a distant observer sees 
one-half of the flux to be radiated. The goal is to show how the apparent angular sizes of 
the masers depend upon the degree of saturation and how, if at all, these sizes depend upon 
the frequency at which they are measured within the spectral line. 

At a specific frequency one-half of the fiux from a spherical maser is emitted within 
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an angle ^1/2(1^) that is found from equation (13) by solving 

f ' I{Re, 9, v) COS 9 sin 9d9^2 f '^'^ ^ 7(i?e, 9, v) cos ^ sin 9d9 (18) 
Jo Jo 

from which it follows that the ratio of the apparent angular size to the geometrical angular 
size is 

^tiJng = sm^9y-,{V) (19) 

The apparent angular size based on the total power is found by performing the additional 
integration over frequency and finding ^1/2 

/•7r/2 /-oo ^ l'0i/2 ^ 

\ I I{R^,9,v) COS 9 sm 9dvd9 = 2 I I I{Re, 9, u) cos 9 sin 9dud9 (20) 
Jo Jo Jo Jo 

to give 

^1/2/^^9 = sin^ ^1/2 (21) 

Results of calculations for ^1/2/^9 as a function of saturation and for Qy^/^g ^ function 
of frequency for representative choices of the saturation are presented in Figures 4 and 5 for 
spherical masers. In addition, the detailed variation of the intensity I{Re, 9, u) with 9 and 
from a representative calculation for a saturated spherical maser is shown in Appendix A. 

The analogous apparent sizes arc given for disks viewed edge-on in Figures 4 and 6. 
Only the angular extent along the plane of the disk varies since the disks here are idealized 
as having zero (or very small and constant) thickness. For disks, 9i/2{i^) is found, using 
equation (14), from 

/ 7(i?e, 9, V) COS 9d9^2 / I{Re, 9, V) cos 9d9 (22) 
Jo Jo 

and then the apparent to geometrical angular extent along the plane of the disk at frequency 
V is 

J^i/2M = sin^i/2(?) (23) 

Similarly, the apparent angular sizes for edge-on, masing disks based on the total power are 
found from 

/ / I{Re,9,u) COS 9dud9 = 2 / / I{R^, 9, u) cos 9dud9 (24) 
Jo Jo Jo Jo 

and 

^1/2/^g = sin^i/2 (25) 
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In Figure 4, we also make comparisons with the results from simpler, approximate calcula- 
tions that are appropriate in limited regimes of saturation (see Goldreich & Keeley 1972). 
In the unsaturated limit, the intensity has the familiar exponential dependence on path 
length. When this form for the intensity is inserted into equations (18) and (22), the appar- 
ent angular sizes can be obtained in a nearly analytic calculation. In the partially saturated 
regime, the reasoning is more involved but still leads to a result that can be evaluated nearly 
analytically. 

In Figures 5 and 6, the normalized spectral line profile of the fiux is shown to demonstrate 
that apparent angular sizes arc relatively insensitive to frequency at frequencies where the 
flux is appreciable. Significant variations in the sizes as a function of frequency do occur. 
However, they are only in the far wings of the line profile. 



The masers in the foregoing Sections 2 and 3 arc isolated. Hence, the seed radiation 
is either the external, isotropic radiation background or the result of spontaneous emission 
within the maser. For some masers, it is likely that a strong continuum source is nearby 
and dominates as the source of the seed radiation. To obtain an indication about how the 
behavior of the spectral line breadths and angular sizes might be different from that of the 
isolated masers, we consider the idealized geometry where the radiation that is incident on 
the maser from the continuum source can be approximated as parallel rays. This idealization 
would be strictly applicable when the continuum source is small and far from the masing 
gas. The rays of maser radiation then remain parallel within the masers and after they 
emerge from the masing region. As a result, there is no interaction between the rays and 
the maser radiation that emerges along a ray is the same as that from a linear maser of the 
same length. The total flux from the sphere or disk is then just the sum of the radiation 
from a series of linear masers with appropriately varying lengths. That is, for the sphere 



4. COMPARISON WITH MASING 
SPHERES AND DISKS NEAR A STRONG 
CONTINUUM SOURCE 




(26) 



where I[L{b),ko,u\ is the solution to 





(27) 



1 + I{s, ko, u) 
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for a linear maser of length L that is analogous to equation (9) but is appropriate when 
the maser radiation propagates in only the one direction because of the strong continuum 
source. Here, L{b) is the path length for a ray that passes through the sphere at impact 
parameter b, and Tq = 2'kqR(, is the unsaturated optical depth parametr for a ray that passes 
through the center of the sphere. For the disk (with the incident rays parallel to the plane of 
the disk), 

F{u,To)^2 I[L{h)XMdh (28) 

^0 

With these expressions for F, the spectral line breadth z/1/2 and the apparent angular sizes 
^1/2 and ^1/2(1^) can be obtained in exactly the same way as was done in Section 3 for the 
F that result from isotropic seed radiation. These are shown in Figures 7, 8 and 9. To 
provide a direct comparison with the linebreadths and the angular sizes for isotropic seed 
radiation in Figures 1 and 4, the linebreadths in Figure 7 and the angular sizes in Figure 8 are 
computed for Jq = 10^^. However, the apparent angular sizes as a function of frequency in 
Figure 9 are obtained for the larger Jq = 10^^ which probably is more realistic when a strong 
continuum source is dominant. As for the case of isotropic seed radiation, the spectral line 
profile in Figure 7 rebroadens in the neighborhood of R/T ~ 1 to the full thermal Doppler 
breadth. Similarly, the variation of the apparent angular size with frequency in Figure 9 is 
unimportant at frequencies within the spectral line profile for which the flux is appreciable. 
The variation of ^1/2 (j^) with frequency is not shown for a sphere. It is similar to that of the 
disk. In contrast, the variation in Figure 8 of the apparent angular size with R/T is quite 
different from that for isotropic seed radiation. The limits for large and small R/T can readily 
be obtained. In the fully saturated regime (large -R/F), /[L(6), /cq, is just proportional to 
the length of the ray with impact paramenter h. This yields 1^1/2/^3 = 0.37 for a sphere 
and i7i/2/^5 = 0.40 for a disk, in good agreement with the angular sizes in Figure 8 at 
logi?/r = 3. In the unsaturated limit, /[L(6), /cq, V\ is proportional to exp(constantxlength). 
For logi?/r — —2, this gives Qi/2/0,g — 0.076 and 0.16, respectively, for the sphere and disk 
— again in good agreement the results in Figure 8. Note that for R/T <^ 1, seed radiation 
due to an isotropic flux and to illumination from one side yield nearly equal Q 1/2/^^9- That 
they should be nearly equal can be seen by inserting the expression for the intensity — with 
its simple, exponential dependence on path length for the unsaturated limit — into equations 
(18) and (22). 



5. DISCUSSION 

The foregoing calculations have been performed for the standard idealization of an 
astrophysical maser — two energy levels and rate equations with "phenomenological" pumping 
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A and decay F rates for the molecular populations as used in equations (1) and (2). This 
standard idealization has been the basis for all previous investigations of the spectral line 
profiles and angular sizes of astrophysical masers that consider deviations from the linear 
maser approximation. The basic equations (6)-(8) that form the starting point are the same 
as those of other investigators (e.g., Bettwieser 1974; Neufeld 1992) as well as for EW. 
The main results are the demonstrations summarized in Figures 1 and 4. (a) The spectral 
line rebroadens to the full Doppler width in the saturated regime for disk and spherical 
masers in essentially the same way as for linear masers and (b) the apparent angular sizes 
of these masers vary only slightly with frequency when the sizes are measured at frequencies 
where the flux is signiflcant ( > 0.1 of the flux at line center). While these conclusions 
were reached in the previous investigation (EW), the solutions were obtained there by using 
approximations to the full equations. That these approximations are of sufficient accuracy 
seems to have been unclear to some readers. In contrast, the calculations here solve the full 
equations that describe the rate equations and radiative transfer for the standard idealization 
for astrophysical masers. In that sense, the solutions are "exact" . Though the calculations 
do depend, in detail, on the intensity of the seed radiation — isotropic background radiation, 
a nearby continuum source, or spontaneous emission, the conclusions (a) and (b) above are 
independent of the type of seed radiation and of its intensity — at least within a plausible 
range of values for the intensity. It should not be surprising that the spectral line rebroadens 
with increasing degree of saturation in a way that is relatively independent of the geometry. 
The maser radiation from a disk or a spherical maser becomes highly beamed when the 
maser is saturated. Highly beamed rays arc nearly parallel, and are thus similar to those of 
linear masers where the rays are exactly parallel. 

All investigations to date of the above issues about the spectral line profiles and apparent 
angular sizes of masers with finite cross sections are based on the same idealized description 
as is used in our calculations. Actual masing regions in astronomy probably deviate at least 
to some degree from these idealized conditions. Velocity gradients, as well as variations in the 
composition and excitation may be present. The observed maser radiation may be a result 
of amplification by spatially separated components. The masing transition may consist of 
hyperfine components (22 GHz water). Finally, A and F are approximations to the results 
obtained by solving a set of rate equations for the velocity-dependent populations of the non- 
masing, as well as the masing, excitation states of the molecule. The rate equations ordinarily 
involve non-masing radiative transitions (usually infrared) and collisional excitations. These 
processes may also cause velocity relaxation without a change in the excitation state of the 
molecule. In view of these complications, the reduction of the rate equations to two-levels 
and phenomenological (A, F) parameters provides a remarkably good description when it has 
been tested (Anderson & Watson 1993). Calculations based on the standard ideahzation thus 
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provide benchmark results. The possibihty cannot be excluded that some of the additional 
complexities that may occur in actual astrophysical masers — which neither we nor others 
have investigated — will lead to deviations from these benchmark results. It seems unlikely, 
however, that the basic conclusion that rebroadening of a spectral line occurs until the full 
Doppler breadth is reached will be altered. 

The observed spectral line breadths of mascr radiation often arc narrower than the 
thermal Doppler breadths for the gas temperatures that are believed to be appropriate for 
the masing gas based on likely pumping schemes or other considerations. Investigators are 
reluctant to accept the interpretation that these masers arc unsaturated because their fluxes 
are relatively constant in time, because the inferred luminosities seem to require a high 
efficiency for the pumping, etc. The natural resolution is that the rate R for stimulated 
emission of these masers satisfies T < R < the rate 7 for velocity relaxation, as recognized 
long ago by Goldreich & Kwan (1974). In this regime, the line breadth remains narrow as 
for unsaturated masers, but the gain (linear with the length of the mascr) and the efficiency 
of the pumping are the same as for a saturated maser. The behavior of the line breadth 
as a function of the degree of saturation has been calculated with velocity relaxation for a 
linear maser (Nedoluha & Watson 1991; also Anderson & Watson 1993) and explicitly shown 
to be the same as without velocity relaxation if 7 is substituted for T. Although velocity 
relaxation has not been incorporated into the calculations of this Paper, its effect on the fine 
breadth is expected to be similar to that for the linear maser. 

We arc grateful to A. Ostrovskii and A. Sobolev for help in finding the best iteration 
algorithms to use. We also acknowledge support from NSF Grant AST99-88104. 

A. Examples of Intermediate Results 

In Figures 10,11, and 12, we show examples of K,{Re,6,u), Jk{Re,w,9') and I{Re,6,u) 
that are computed at the surface of a masing sphere that is radiatively saturated. In the 
calculations, quadratic interpolation is used for k and near the symmetry points where 6 
or 9' is equal to or 90 degrees, and at u and w = 0. 

B. The Geometry for the Calculation of Isolated Disk and Spherical Masers 

To perform the computation with a disk or spherical maser for the intensity /(r, 9, v) at 
radius 1^ propagating in the direction k (which makes an angle 9 with 1^) and with frequency 
I/, equation (6) is integrated along the paths of the rays such as the ray shown in Figure 
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13. For example, at the integration step that is a distance s along this path with impact 
parameter b shown in Figure 13, k from the previous iteration is evaluated at the location 
"r*" for the angle 6" to compute the right hand side of equation (6). After the intensities at 
for rays at a sufficient number of angles have been computed, the integral in equation (8) 
is performed to find Jk at the location for the molecular velocity of magnitude w that is 
in the direction 6' shown in Figure 13. In turn, these values of Jk{r,w,6') are inserted into 
equation (10) for a spherical maser or equation (12) for a disk maser to obtain a new value 
for K{r, 9, v). In the computation for a spherical maser, ■uf ordinarily is not in the same plane 
as V and k. The V/c plane differs from the 1^1$ plane by the angle of rotation (f) about 'r' 
as the axis of rotation. 



C. Additional Information on the Numerical Methods 

For the reader who might be interested in repeating our numerical calculations, here 
are a few details. In order to place more integration points near the center of the sphere 
(or disk), where as discussed by, e.g. Goldreich and Keeley (1972), there is a transition 
region between saturated and unsaturated, we distributed the points at radii proportional to 
r^. Recognizing that the radiation is beamed, we distributed the angular points uniformly 
over the range < < ^^max with ^^^ax = 180°/(1 + 0.2/s:or). Thus, 9 covers the entire 
angular range at the center, and decreases toward the perimeter and with saturation since 
hqt increases with increasing saturation. The coefficient 0.2 in the expression for ^max was 
chosen by trial and error. We also experimented with other forms for ^max and for the r 
dependence. We examined the choices which we adopted for the distributions of points in 
angle and in radius by plotting various distributions of the computed quantities after the 
fact. 

The equations were solved by cross iteration between the equation for Jk and the equa- 
tion for K using the linear superposition Jk — fJl + (1 — f)Jl to obtain Jk for the next 
iteration step from the results of the most recent iteration step (J^) and the next most re- 
cent iteration ( J^) step. We used / = 0.4 for the first 15 iterations, / = 0.2 for the next 15, 
and / = 0.1 for the subsequent iterations which typically number about 100 for high satura- 
tion (with fewer intcrations for lower saturations) to achieve the convergence criterion that 
the maximum change between iteration steps in Jk and k at any integration point be less 
than one percent. Further tightening of the convergence criterion led to negligible change in 
the results. Computations are performed first for a completely unsaturated case (small noRe) 
where Jk ^constant is taken as the initial value for the iterations. For subsequent, larger 
values of kqR^ (where saturation may be important), the converged Jk from the computation 
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for the previous choice for KoRe are taken as the initial values after rescaling in r. 

The calculation did not require a supercomputer, even for the case of the spherical 
geometry which is the most difficult. All computations were performed on a PC with an 
AMD XP1900 microprocessor running under MS Windows. When the code is compiled using 
Compaq Fortran, the numerical computations for single curve in Figures such as Figure 
1 requires about 36 hours for a sphere. When the code is compiled by Microsoft Visual 
Studio. Net optimizing C++, the analogous running times were typically 30% less. 

REFERENCES 

Anderson, N., & Watson, W.D. 1993 ApJ, 407, 620 
Bettweiser, E.V. 1974 A&A, 50, 231 
Bettweiser, E.V. & Kegel, W.H. 1976 A&A, 37, 291 
Emmering, R.T., & Watson, W.D. 1994 ApJ, 424, 991 
Goldreich, R, & Keeley D.A. 1972 ApJ, 174, 517 
Goldreich, R, & Kwan, J.Y. 1974 ApJ, 190, 27 
Gwinn, C.R. 1994, ApJ, 429, 253 
Litvak, M. 1971 ApJ, 170, 71 

MoscadeUi, L., Menten, K. M., Walmsley, C. M., & Reid, M. J. 2003 ApJ, 583, 776 

Nedoluha, G.E., & Watson, W.D. 1991 ApJ, 367, L63 

Neufeld, D.N. 1992, ApJ, 393, L37 

Wallin, B.K., & Watson, W.D. 1997 ApJ, 476, 685 

Watson, W.D., Sarma, A.R, & Singleton, M.S. 2002 ApJ, 570, L37 



This preprint was prepared with the AAS IM^jX macros v5.0. 



-18- 




0.2 - 



' ^ ^ ^ ^ ' 

-2-10123 

iog(R/r) 

Fig. 1. — Comparison of the spectral line breadth (FWHM) divided by the thermal Doppler 
breadth (FWHM) for the flux of observed maser radiation from linear, thin (edge-on) disk 
and spherical masers as a function of the log of the degree of saturation R/T (with Iq — 10"^). 
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Fig. 2. — (upper panel) The quantity b which measures the deviation of the profile of the 
radiation from the Gaussian that best fits the profile, as a function of the log of the degree of 
saturation, for the hnear, thin (edge-on) disk and spherical masers (with /q = 10~^). (lower 
panel) The kurtosis K for the same profiles. 
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Fig. 3. — The dependence of the spectral hnebreadth of the observed maser radiation on Iq 
for Iq = 10-^ 10-^,and 10"^ The spectral hne breadth (FWHM) divided by the thermal 
Doppler breadth (FWHM) for the flux of observed maser radiation from linear, thin disk, 
and spherical masers is given a function of the log of the degree of saturation. 
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Fig. 4. — The log of the ratio Qi/2/f2g of the apparent angular size to the geometrical 
angular size for thin (edge-on) disk and spherical masers as measured by the total power, as 
a function of the log of the degree of saturation (with /q = 10~^) . Apparent angular sizes 
based on approximations that are appropriate in the unsaturated and partially saturated 
regimes are given by the straight lines (see text). 
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Fig. 5. — The ratio VL^j^/^g of the apparent angular size to the geometrical angular size for 
a spherical maser as measured by the flux at frequency v. The ratio is shown for values of 
the degree of saturation R/V {— 4, 10, 100) as indicated (with /q — 10~^). The normalized 
flux (normalized by the flux at z/ = 0) also is shown as a function of frequency at the same 
values of R/T. 



-23- 




frequency 



Fig. 6. — Similar to Figure 5 except the maser is in a thin disk. 
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Fig. 7. — Same as Figure 1 except that a strong, external continuum source on one side (the 
far side) provides the seed radiation for these masers. The spectral line breadth (FWHM) di- 
vided by the thermal Doppler breadth (FWHM) is shown for linear (dots), thin disk (dashes), 
and spherical (solid line) masers. For comparison purposes, the spectral line breadths in Fig- 
ure 1 for an isotropic background are shown as well. The background continuum intensity 
for all of the computations in this Figure is /q = 10~^. 
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Fig. 8. 



Similar to Figure 4 except that a strong, external continuum source on one (the far 
side) provides the seed radiation for these masers. The ratio of the apparent angular 

size to the geometrical angular size for the thin disk and spherical maser as measured by the 
total power is shown. The ratio is given as a function of the log of the degree of saturation 
R/r. For comparison purposes, the analogous ratios fli/2/^g from Figure 4 for thin disk and 
spherical masers with an isotropic background are shown as well. The background continuum 
intensity for all of the computations in this Figure is Iq = 10~^. 
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Fig. 9. — Similar to Figure 5 except that a strong, external continuum source on the far 
side provides the seed radiation for these masers. The ratio ^^j^j^g of the apparent angular 
size to the geometrical angular size is shown for a thin disk maser as measured by the 
flux at frequency v. The ratio is given for several values of the degree of saturation RjV 
{— 1, 10, 100) as indicated (with /q — 10~^). The normalized flux (normalized by the flux at 
= 0) also is shown as a function of frequency at the same values of RjV. 
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Fig. 10. — An example of the function Jk{Re, w, 9') in the computations for a spherical maser 
with seed radiation due to an isotropic background (Jq = 10~^). In the computations for 
this example, R/T — 52. 
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Fig. 11. — An example of the function li{Re, 6, v) in the computations for a spherical maser 
with seed radiation due to an isotropic background (/q = 10~^). In the computations for 
this example, R/V = 52. 
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Fig. 12. — An example of the intensity I{Re,6,u) at the surface in the computations for 
a spherical maser with seed radiation due to an isotropic background (Jq = 10~^). In the 
computations for this example, R/T = 52. Note that numerical values on the intensity axis 
should be multiplied by 10^. 




Fig. 13. — Diagram of the geometry. 



